Glauber slow dynamics of the magnetization in a molecular Ising chain 
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The slow dynamics (10~^ s - 10* s) of the magnetization in the paramagnetic phase, predicted 
by Glauber for Id Ising ferromagnets, has been observed with ac susceptibility and SQUID magne- 
tometry measurements in a molecular chain comprising alternating Co^"^ spins and organic radical 
spins strongly antiferromagnetically coupled. An Arrhenius behavior with activation energy A=152 
K has been observed for ten decades of relaxation time and found to be consistent with the Glauber 
model. We have extended this model to take into account the ferrimagnetic nature of the chain as 
well as its helicoidal structure. 

PACS numbers: 75.1Q.H, 76.90.+d, 75.40. Gb 



Slow magnetic relaxation in low dimensional systems 
without long range order is a very interesting subject 
both from the fundamental and the applicative point of 
view This phenomenon has been clearly observed 
and theoretically explained in molecular magnetic clus- 
ters, which show also hysteretic effects of molecular ori- 
gin 1^. These clusters can be considered as magnetic 
quantum dots, i.e. zero-dimensional compounds. In one- 
dimensional (Id) systems, a macroscopically slow relax- 
ation was predicted many years ago by Glauber |^ for a 
ferromagnetic Ising chain. In the last three decades many 
quasi-ld magnetic compounds were synthesized and char- 
acterized and new features were observed, such as non- 
linear excitations Haldane gap § and spin-Peierls 
phase j^. On the contrary, up to now, Glauber's predic- 
tion has not been experimentally confirmed because two 
strict conditions must be met: i) a very strong easy axis 
anisotropy, ii) a very high ratio J/ J', between the intra- 
chain and the interchain exchange interactions, leading 
to a very low critical temperature to three-dimensional 
(3d) order. 

We have recently reported preliminary data which sug- 
gest that these conditions can be fulfilled in a real quasi- 
ld ferrimagnetic compound, [Co(hfac)2NITPhOMe] 
(CoPhOMe in the following) ||^. It has a trigonal sym- 
metry with ternary helices [||, along the c axis, formed 
by Co(hfac)2 moieties and by the NITPhOMe organic 
radicals, as outlined in Fig. 1. The primitive cell con- 
tains three Co^+ metal ions alternating with three rad- 
icals. Each Co^+ is geometrically related to the other 
two through a 120° rotation around the c axis and has 
a local axis of easy anisotropy, z, which makes an angle 
9 ~ 50° with c, as observed in monomeric complex of 
formula Co (hfac) 2 (NITPhOMe) 2 [§, where the Cobalt 
ion coordination is practically the same as in CoPhOMe. 
The strong Id character is a consequence of the strong 
metal-radical magnetic exchange along the chain and of 



the absence of chemical bond between different chains. 
The temperature dependence of xT gives a clear indi- 
cation of Id behavior, with antiferromagnetic coupling 
between 00^+ and NITPhOMe radical. Since the effec- 
tive 5 = 1/2 spin of the Co^"*" and the S = 1/2 spin of 
the radical are characterized by different g values {gco 
and gji in the following), the chains behave as Id ferri- 
magnets, due to the non-compensation of the magnetic 
moments. An approximate quantitative analysis using 
the alternating spin Ising model suggested that the in- 
trachain exchange interaction is of the order of 200 K 
in absolute value. The ac magnetic susceptibility shows 
that the relaxation time of the magnetization below 20 
K follows a thermally activated behavior. Below 4 K a 
stepped hysteresis was observed with the field parallel to 
the trigonal axis, while no hysteresis was observed with 
the field perpendicular to c. 

We wish to work out here an extension of the Glauber 
dynamics to a helical Ising ferrimagnet, showing how it 
can unambiguously justify the experimental features re- 
called in more detail here below. In Fig. 2 the isothermal 
magnetization as a function of the field H is reported at 
T = 2 K for two different orientations of the field. When 
H is applied along the c axis, a hysteresis cycle appears 
below 5 K, while the hysteresis is absent down to 2 K for 
fields applied perpendicularly to c. The unusual shape 
of the hysteresis (characterized by well defined steps at 
H = ±4 kOe) is discussed later on. The magnetization, 
measured by cooling in zero field and then applying a 
field of 100 Oe, shows the onset of irreversibility around 
6 K, if measured with the field parallel to the c axis, while 
a paramagnetic behavior is observed down to 1.8 K by 
applying the field in the trigonal plane. The blocking of 
the magnetization at 6 K is not due to a phase transition 
to 3d order but has a dynamical nature, as confirmed 
by ac magnetic susceptibility measurements performed 
in the range 0.2 Hz - 95 kHz, which show a strong fre- 
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quency dependence below 20 K: the in-phase suscepti- 
bility x' (not shown) goes through a maximum and the 
out-of-phase component %" becomes different from zero 
and goes through a maximum (see Fig. 3) , showing that 
the relaxation in this temperature range becomes slow in 
the time scale of the experiment (10^^ s - 1 s). 

The observed relaxation time r follows an exponential 
law due to a thermally activated mechanism (see Fig. 4) 
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The best fit gives to = 4(1) • 10"" s, A/kg = 152(1) 
K. The value of tq is significantly smaller than that ob- 
served in superparamagnetic clusters |^ . It is worth 
noticing that x' goes to zero for T below the peak in x" 
and indeed x' 8.nd x" almost describe a semicircle in the 
Cole-Cole plot (see the inset of Fig. 3) ||ll|]. The ac suscep- 
tibility recorded with the field oscillating in the trigonal 
plane is essentially paramagnetic except for a small frac- 
tion, which we attribute to the non perfect alignment of 
the bunch of oriented crystals used for these measure- 
ments. For T < 5 K the relaxation was also studied by 
monitoring the decay of the longitudinal magnetization 
which was found to be exponential, as a function of time, 
if weak fields (< 1 kOe) were applied. The characteristic 
time T is in agreement with the law (1) and the parame- 
ters reported above. The Arrhenius behavior is therefore 
obeyed in the wide temperature range 4.5 K - 15 K and 
over ten decades in time (see Fig. 4). In a real system in- 
terchain interactions are always present, however the 3d 
ordering cannot be detected if the critical temperature is 
well below the magnetization freezing temperature. This 
hypothesis is supported by the lack of anomalies in the 
specific heat |12 and in preliminary neutron diffraction 



spectra |13 



The dynamic behavior is therefore strongly reminiscent 
of that observed in molecular clusters |^ [|o| , where the 
energy barrier which hampers the reversal of the mag- 
netization is originated by the magnetic anisotropy. In 
the present case, despite the large value of the barrier, 
the plane orthogonal to the trigonal axis is not a hard 
plane, as the magnetization rapidly approaches satura- 
tion (see Fig. 2). This anomaly appears to be related 
to the helicoidal structure of the chain, as demostrated 
by the following extension of the Glauber theory to a 
helicoidal ferrimagnetic chain. 

According to the Glauber model the j-th spin aj 
is assumed to take the integer values ±1, corresponding 
to two possible projections along its local z axis. The 
aj coordinates depend stochastically on time. If we con- 
sider an isolated spin, the probability per unit time to 
reverse is given by |ao: where ap is a free parameter 
of the model. In the case of N coupled spins we need 
another phenomenological parameter, 7, which accounts 
for the tendency of the j-th spin to align parallel to its 



nearest neighbors: i.e., the transition probability from 
the state aj to the state —a-j is assumed to be 
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^ao[l - ^Wj{aj+i + cTj-i)] (2) 



Imposing the detailed balance it is possible to obtain 
a correspondence between 7 and the exchange constant 
J of the Ising model Hamiltonian. When the system 
has reached equilibrium at temperature T, the ratio 
of the transition probabilities Wa-j^-a-j and w-^j^oj 
must equal the ratio of the Maxwell-Boltzmann prob- 
abilities associated with the two equilibrium configu- 
rations (fJi, -CTj, (Tat) and (di, ctj, ctat) 
This yields 7 = tanh( J/2A:bT). 

In the Glauber model the direction of each local z axis 
is determined by the applied field while in our case we 
will consider as z axis of each Co^+ the local axis of easy 
anisotropy. So inside the same cell we will have three dif- 
ferent z axes for the Co^^ ions (z2, Z4, zq) and three for 
the radicals (zi, Z3, Z5, which are given by the vectorial 
sum of the directions of magnetic anisotropy of the two 
Co^^' nearest neighbours). The extended Ising hamilto- 
nian of this system of spins can be written: 
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where is the unitary vector corresponding to the Zj 
direction. In the coordinate (T„j , n represents the cell 
index and j the magnetic center inside the n-th cell. 
Since we are interested in the ac response of the sam- 
ple to an applied oscillating field, we can focus on the 
stationary solution of the problem. So the system of 6A^ 
coupled equations, describing the Markov relaxation pro- 
cess of the spins chain, reduces to six coupled differen- 
tial equations for the expectation values Zj{t) = {un,j{t)) 
(j = l,---,6). In fact the translational invariance al- 
lows to neglect the n index. The equation describing the 
evolution of the system can be written in matricial form 
as 
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The k vector depends on the orientation of the oscillating 
magnetic field with respect to the local axis of each spin 
and its elements will be specified later on. The symmet- 



ric matrix A has elements Aj 



-1, A, 



Ai.g = Aq^i = 7/2 and zero otherwise. 
Denoting by U the orthogonal matrix which diagonalizes 
A, D — U^AU, we perform the transformations ( — U^z, 
K — U^k, so that the system (4) is decoupled as 
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where the diagonal matrix D has elements Dn = —1 +7, 

D22 = -1 - 7, D33 = Dm = -1 + 7/2, D55 = Dee = 
— 1 — 7/2. Each equation of the decoupled system (5) can 
be integrated, in the stationary limit, to give 



(6) 



with tti = (1 - 7)q;o; 02 = (1 + 7)ao; as = 04 = 
(1 — -f/2)ao; a5 — ae — {1 + j/2)aQ. Since I7I < 1, we 
have that 03, a^, as, ag cannot vanish. On the contrary, 
for T — > 0, ai and a2 can vanish for J > and J < 
respectively. In our case J < and the slow relaxation 
(exponentially diverging relaxation time) is expected for 
\J\/kBT > 1 due to the mode (2 
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while Ti ~ l/tto- The law (7) is in agreement with the 
fit (1) where A = —J and tq = i^^. 

When the magnetic field is applied along the c axis, its 
projection along each local axis e^ is trivially H{t)cos9 
and, while K3|| = K4|| = K5|| = Kgy = 0, one has 
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where rj = tanh ( 
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For Kill the difference between 



the Lande factors is substituted by the addition. The 
complex susceptibility can be obtained through the pro- 
jection of the contribution coming from each spin: 



X (w) = S'Vb cos^ 9 
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where, since in the experimental conditions (frequency 
and temperature) <C cki, the first term is negligible. 
Hence, by varying the frequency, x" as a function of %' is 
expected to describe a semicircle, in agreement with the 
experimental data reported in Fig. 3. If the oscillating 
magnetic field is applied along a generic perpendicular 
direction with respect to the c axis, after a few steps we 
find = K2_L — and consequently the slow mode C2 
is not activated but only a linear combination of (^3, (^4, 
Cs , Ce is. These results explain the experimental evidence 
for slow relaxation for magnetic fields applied parallel to 
the c axis while, for magnetic fields perpendicular to c, 
the relaxation time r remains of the order of — . 

The hysteresis curves reported in Fig. 2 are character- 
ized by well-defined steps, the more evident observed at 
H = ±4 kOe. The steps observed in the first and third 
quadrants, i.e. when the magnetization does not change 
sign, have a static nature, as confirmed by the fact that 
they are fully reversible (see inset of Fig. 2), do not 



depend on the sweeping rate of the field, and are observ- 
able above the blocking temperature, when the hysteresis 
disappears. Therefore, their origin seems to be differ- 
ent from the resonant tunneling mechanism observed in 
molecular magnetic clusters Q Similar steps of 

static nature were observed in linear antiferromagnetic 
chains with alternating magnetic moments when the Zee- 
man energy is comparable to the exchange interaction 
]l6[ . In the present case we expect the exchange interac- 
tion to be much larger than the Zeeman energy in a field 
of 4 kOe. However the helical structure and the fact that 
the local z axes are not pointing along the axis of the helix 
[9 ~ 50°) can account for this anomaly as well as for the 
similar saturation values of the two components of the 
magnetization (Fig. 2) despite the completely different 
dynamic behavior. This hypothesis has been confirmed 
by classical transfer matrix calculations and by diagonal- 
ization of a finite quantum ring with six spins, where the 
trigonal symmetry of the helical structure has been taken 
into account. A detailed description of this study will be 
the subject of a future publication. 

In conclusion, we have shown that CoPhOMe is a 
quasi-ld ferrimagnetic compound which at low tempera- 
ture shows slow dynamics of the macroscopic magnetiza- 
tion of the type predicted many years ago by Glauber and 
never observed before. The presence of magnetic hystere- 
sis in the absence of 3d magnetic order opens the exciting 
possibility of storing information in an array of magnetic 
centers with a dramatic reduction in the dimensions (two 
of the three) of the magnetic memory unit. From a more 
fundamental point of view, we have extended the the- 
ory developed by Glauber for the simple model of an 
Ising ferromagnetic chain to the case of helicoidal geom- 
etry, and the anomalous dynamic behavior of CoPhOMe 
turned out to be completely understandable in terms of 
it. 
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FIG. 1. Schematic view of the helicoidal structure of 
CoPhOMe. The laboratory reference frame corresponding to 

the crystallographic axes a* ,b, c and the local reference frame 
(x, y, z) are shown. The angle between z and the c axis is 
the angle mentioned in the text. 



FIG. 2. Hysteresis loops recorded on a single crystal of 
CoPhOMe by applying the field in the a*h plane (solid line) 
at T = 2 K and along the c axis at 2.0 K (black squares), 3.5 
K (open circles) and 4.5 K (open triangles). In the inset: at 
T = 2 K the field was cycled three times around the step at 
4 kOe. The curves taken with increasing field (solid circles) 
and decreasing field (solid line) are superimposable, denoting 
the reversible character of the step. 

FIG. 3. Temperature dependence of the imaginary compo- 
nent of the ac magnetic susceptibility measured in zero static 
field and with the dynamic field oscillating parallel to the c 
a:xis for nine selected frequencies in the range 0.18 Hz - 95 
kHz. In the inset x" vs x' (i-e, the Cole-Cole plot) measured 
at 10 K. 

FIG. 4. Temperature dependence of the relaxation time. 
The data represented by circles have been obtained from the 
curves of Fig. 3 by assuming that, at the temperature of the 
maximum in each curve, r = 1/{2ttu), where v is the fre- 
quency of the oscillating field. The triangles represent the 
relajcation time obtained from the decay of the magnetization 
after having applied and removed a weak field. The solid line 
represents the fit of the x" a-c data using Eq. (1). The best 
fit parameters are A = 152(±1) K and to = 4(±1)10"" s. 
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